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Abstract. The main aim of this paper is to discuss recent results on the adaptation 
of the Fokas-Gel'fand procedure for constructing soliton surfaces in Lie algebras, which 
was originally derived for PDEs [Grundland, Post 2011], to the case of integrable ODEs 
admitting Lax representations. We give explicit forms of the g-valued immersion 
functions based on conformal symmetries involving the spectral parameter, a gauge 
transformation of the wave function and generalized symmetries of the linear spectral 
problem. The procedure is applied to a symmetry reduction of the static 4 -field 
equations leading to the Jacobian elliptic equation. As examples, we obtain diverse 
types of surfaces for different choices of Jacobian elliptic functions for a range of values 
of parameters. 
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1. Introduction 

The methodological approach assumed in this work is based on the Fokas-Gel'fand 
formula for immersion of 2D-soliton surfaces associated with integrable models. This 
topic has been extensively developed by several authors (see e.g. [H El El HS1 IE]) 
and also by the authors [8] for integrable partial differential equations (PDEs). The 
results obtained proved to be so fruitful from the point of view of constructing 2D- 
surfaces immersed in Lie algebras that it seemed worthwhile to adapt this method and 
check its effectiveness for the case when an integrable ODE can be written in the Lax 
representation. The motivation for such an analysis is that it allows for the study of 
PDE surfaces using ODE surfaces as approximations. The construction of such soliton 
surfaces using the Fokas-Gel'fand approach is presented in [H] and its application to the 
symmetry reduction of the static 4 -field equations is the goal of this work. 

2. Soliton surfaces associated with integrable ODEs 

Let us consider an ODE (x stands for the independent variable and u for the unknown 
function) 

A[u] = A(x,u, u x , u xx , . . .) = (1) 

which admits a Lax pair with potential matrices L(\, [u]), M(A, [u]), taking values in a 
Lie algebra q, which satisfy 

D X M + [M, L] = 0, whenever A[u] = 0, (2) 

where L(A, [it]), M(A, [u]) are rational functions of a spectral parameter A taking values 
in either C or R. 

In what follows, we make use of jet space and the prolongation structure of vector 
fields as presented in the book by P. J. Olver [T3]. For derivatives of u we use the 
standard notation u x and uj for the first and Jth derivatives of it with respect to x, 
respectively. Functions depending on the independent variable x, dependent variable u 
and its derivatives are denoted 

f[u] = f(x,u,u X} u XX} . . .) 

and the total derivative in the direction of x takes the form 

d d 

D x = d x + u x — + u xx - h.... (3) 

au au x 

This Lax pair equation ([2| can be regarded as the compatibility conditions of a 
linear spectral problem (LSP) for the wave function $(A, y, [u]) taking values in the Lie 
group G 

D x $(\,y, [it]) = L(X, [u))®(\,y, [«]), D v $(\,y, [u]) = M(A, [«])$(A,j/, [it]). (4) 
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Note that the wave function <3> depends on an auxiliary variable y in the LSP Q and so 
the compatibility condition for Q coincides with ^ since u y = 0. The total derivative 
in the direction y is given by 

and consequently, we obtain 

D y L = D y M = 0. (5) 

It was shown in |6J, that for any g-valued functions A(X,y, [u]) and B(X,y,[u}) 
which satisfy 

D y A - D X B + [A, M] + [L, B]=Q (6) 

there exists a g-valued function F with tangent vectors given by 

D X F = D y F = (7) 

Whenever the matrices A and B are linearly independent, F is an immersion function 
for a 2D surface in the Lie algebra q. It was proved in [6j [71 [8] that there exist three 
linearly independent terms satisfying ^ given by 

A — a(X)—L + D X S + [S, L] +prv Q L e g, (8) 

B = a(\)—M + D y S + [S, M] + prv Q M G g, (9) 

where a(A) is an arbitary scalar function of A, S = S(X,y, [u]) is an arbitrary g-valued 
function, vq is a generalized symmetry of ([TJ and its prolongation is given by 

v Q = Qlu} — , prv Q = v Q + Dj(Q) — . (10) 

Here vq is assumed to be a generalized symmetry of a nondegenerate ODE ([!]), i.e. 

Ptvq(A[u]) = 0, whenever A[u] = 0. (11) 

Further, it was shown in [8] that the g-valued function F can be integrated (up to an 
additive g-valued constant as) 

F = a(A)$" 1 ^ + <&- l S& + Q^prvq® e 0, (12) 
o X 

as long as vq is a generalized symmetry of the LSP Q and the ODE in Lax 
representation (fij) in the sense that the following equations hold 

prv Q (D X M + [M, L]) = 0, whenever D X M + [M, L] = 0, (13) 
prv Q (D x $ - L$) = 0, whenever D x $ - L$ = 0, (14) 
prv Q (D y ® - M$) = 0, whenever - M$ = 0. (15) 

The three terms in (|6| correspond to a conformal transformation in the spectral 
parameter (known as the Sym-Tafel formula for immersion [121 US]), a gauge symmetry 
of the LSP (|4]) and generalized symmetries of the ODE ([!]) and the LSP Q. The 
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integrated form (12) define a mapping F : IR 2 — > q which is called the Fokas-Gel'fand 
formula for immersion in a Lie algebra q. In what follows, we will refer to it as such. 

As first indicated in [6] and extended in [TO], there are many more choices of A 
and B which satisfy (j6|. In fact, as proven in [10], any q- valued function on jet space 
can be transformed into a symmetry of the Lax equation ^ when considered as an 
autonomous system of PDEs in the variables L and M. For example, this system is 
invariant under the following point transformations: translation in the x direction and 
conformal transformations in the spectral parameter A corresponding to terms associated 
with scalar constants ot\ and a 2 in ( Jl6| ) and (17). The symmetry of the equation with 
respect to an arbitrary gauge S(X,y,[u\) G Q for the wave functions $ of the LSP 
Q is given by the term associated with a^. The expressions corresponding to and 
a 5 are related to the invariance of Q under dilations (i.e. x — > e^x, L — > e~^L and 
y — > e^y^M — > e~ M M, /j 6 1). This leads to the following extension of the matrices 
A,B EQ given by g and g to 
dL 

A = axD x L + "27^- + a 3 (D x S + [S, L\) + a 4 D x (xL) + a 6 prv Q L, (16) 
dM 

B = a 1 D x M + a 2 ^— + a 3 (D y S + [S, M)) + a 4 xD x M + a 5 D y (yM) + a 6 prv Q M. (17) 

OA 

The case where «i = 0:3 = «4 = a$ = a§ = 0, «2 = a(X) corresponds to the 
Sym-Tafel formula for immersion [151 US] > which is given by 

F 5T = a(A)$- 1 ^$G S (18) 
with tangent vectors 

D X F ST = a(A)$- 1 ^L$, D y F ST = a(A)$- 1 ^M$. (19) 



If the tangent vectors (19) are linearly independent, then the function F is an 
immersion of a 2D-surface in the Lie algebra g. 

The case where «i = «2 = «4 = «5 = «6 = was studied in [U [7J E] and the 
corresponding surfaces can be integrated explicitly as 

F s = ^- 1 S(X,y,[u])^eQ (20) 

with tangent vectors 

D X F S = $- 1 (D x S+[S,L})$, D y F s = $- 1 (D y S+[S,M})$. (21) 

For F s to be an immersion, we require the linear independence of the tangent vector 
fields. 

Here, we show only the proofs for the case involving translation in the x direction 
and the terms involving dilations. In the first case, (when ct 2 — «3 — «4 — «5 — «6 — 0), 
the tangent vectors are given by 

D X F = <5>- 1 (D X L)<!> 1 A = D X L } (22) 
D y F = <b-\D x M)§, B = D X M. (23) 
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The matrices A and B satisfy the condition ^ 

D y (D x L) - D X (D X M) + [D X L, M] + [L, D X M) 

= -D x (D x M+[M,L}) = 0, 

whenever the Lax equation ^ holds. Thus, this proves that there exists a g-valued 
function F with tangent vectors given by (22) and (23). Furthermore, the immersion 
function can be integrated and is given by 

F = ^D x ^ = (24) 

whenever the wave function $ is a solution of the LSP Q. It is straightforward to check 
that 

D X F= -$- 1 (D x $)$- 1 L$ + $- 1 D x (L$) 



*-\D x L)$, 



and 



DyF 



= -$ _1 [L,M]$ 
= $- l (D x M)$, 

whenever (f2l) holds. Here we have used D y L = 0. Hence, we have proved that the 
immersion function F can be integrated as (24). 

In the case of dilation symmetry (where ai = «2 = «3 = «6 — 0) the tangent 
vectors have the form 



with 



D X F = D V F = 



A = a 4 D x (xL), B = a 4 x(D x M) + a 5 D y (yM). 



(25) 



(26) 



It is a straightforward computation that matrices A and B (26) satisfy the condition 
(JgJ) whenever the Lax equation holds ([2]). From (25), we can integrate and find the 
immersion function 

F = a A x^~ x L^ + a 5 s/$ -1 M$. (27) 

So we have, 

D X F = a^- 1 (L - xL 2 + xD x L + xL 2 ) $ + a^ 1 (-LM + D X M + ML) $ 
= a 4 &~ 1 D x (xL)&, 

and 

D y F = a A x (-ML + D y L + ML) $ + oc b ^~ l [M - yM 2 + yD y M + yM 2 ) $ 
= a 4 x<$>- l D x M<$> + a^DyiyM)®, 
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whenever the wave function <3> satisfies the LSP Q and the potential matrices satisfy 
the Lax equation ([2]). This shows that the immersion function F can be integrated as 
(27) and that the vector fields 

v Qa = D x (xL)—+x(D x M) — , (28) 

and 

d 

vn = M— (29) 

are generalized symmetries of the Lax equation ([2]). 

Finally, we demonstrate the proof of the integrated form of the surface for the final 
term (when ati — . . . — 015 = 0), as shown in a similar way in [8j. Suppose that vq is a 
generalized symmetry of an ODE written in Lax representation ([2]). Then, the matrices 

A = prv Q (L), B = prv Q (M), (30) 

satisfy condition (J6| 

D y A - D X B + [A, M] + [L, B] = 

= D y {prv Q L) - D x (prv Q M) + [prv Q L,M] + [L,prv Q M] 

= prv Q (D y L - D X M + [M, L]) . (31) 

if and only if vq is a generalized symmetry of ([2]). Here, we have used the fact that a 
generalized vector field in evolutionary form commutes with total derivatives [H] 

\pr0 Q ,D x ]=O, [prv Q ,D y }=0, (32) 

and that the potential matrices L and M do not explicitly depend on y (ie. @ holds). 
Thus, there exists a g-valued immersion function F for a surface with tangent vectors 

D X F = $-\prv Q )$, D y § = $- l (prv Q M)$, (33) 



as long as the matrices A and B (30) are linearly independent. Furthermore, the 
immersion F can be integrated explicitly as 

F = <&- l prv Q §, (34) 
if and only if the vector field vq is a generalized symmetry of the LSP Q. That is, 

D X F = -$- 1 L(prv Q $) + $- 1 D x (prvQ$) 
= - ^Liprvq®) + ^- 1 prv Q {D x ^) 
= - $ _1 prwQ(L$) + $ _1 (pri?QL)$ + ^ 1 prvg(D x ^) 
= $ _1 (pr-UQL)$, whenever prvQ (D x & — L$) = 0, 

and the proof is similar for D y F. Thus, the immersion function F can be integrated as 
(34) if and only if the vector field vq in evolutionary form is a generalized symmetry of 
the LSP Q and ODE in Lax representation ([2]). 

Hence, for the system composed of equations (|2]) and Q, the corresponding formula 
for immersion (12) becomes (up to an additional q- valued constant) 

F = (^onD x + a 2 a(A) + a 3 S + a 4 xL + a 5 yM + prv Q ^j $, (35) 
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where vq is a generalized symmetry of both equations ^ and Q. 

In the next section, we consider second-order autonomous ODEs and their Lax 
pairs. The wave function for the associated LSP is given explicitly and can be used for 
the purpose of constructing soliton surfaces. 



3. Second-order autonomous equations 



Let us consider a second-order, autonomous differential equation given by 



u , 



f(u), f'(u) 



d_ 
du 



for some function f'(u). Equation (36) admits the first integral 



u , 



and its solutions are known and satisfy 
du 



6{X - Xq 



Xq E 



(36) 
(37) 
(38) 



Note that, in the case when f(u 



,-1/2 



R(u, y/Pju)) is a rational function of its 



arguments and P(u) is a polynomial of degree 3 or 4, the function u which solves 



(36) is the inverse of an elliptic integral [2]. Note also that the constant of integration, 



in the first integral (37), can be absorbed since the function / is an arbitrary function 
of u. 



The ODE ( 36 ) admits the following Lax pair 



D x $ 

where the potential matrices 





f'(u) 
u+X 







M 



Ux 

u + X 



/M-g(A) 

u+X 

-u x 



(39) 



take values in the Lie algebra si (2, M.) and are rational functions of the spectral parameter 
whenever g(X) is a rational function of A. Note that det(M) = —g(X) and the choice 

<?(A) = /(-A) 

makes the potential matrices L and M polynomial in u whenever f(u) is a polynomial 
in u. In what follows, we call g(X) the discriminant. 

The goal is to construct surfaces in the Lie algebra q by the Fokas-Gel'fand 



procedure for the general form of the ODE (36). For this purpose, we solve the LSP 



Q and find explicitly the most general form of the wave function $ = G SL{2 
with components [9] 

M2+ + 02-), $22 : — =\<.>2- 



$21 



12 



2v^ 
1 



V2' 



2V9 



(40) 
(41) 
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where 



-'lit 



2± = Vu + A*±, 



± 



exp 



±VS(A) 



2/ 



2(m + A)y_ ' 
The generalized vector field in evolutionary form 
9 



du 



is a generalized symmetry of the ODE (36) if and only if 

= 0, 



prv Q u 



whenever u xx — -f'(u) = 



holds. The determining equation for Q is 



1 



D X Q — -f"(u)Q = 0, whenever u xx — -f'(u 



0. 



(42) 

(43) 
(44) 

(45) 
(46) 
(47) 



2 v ' T 
Here f\u) and f"{u) are the first and second derivatives of / with respect to u. It is 
straightforward to verify that Q = u x is a solution of the determining equations (47). 
Therefore, for an arbitrary function f(u), the vector field v Ux is a symmetry of ([2j and 
Q, since the the prolongation of acts as a total derivative on functions which do 
not depend explicitly on x, as is the case for both (|2j) and Q. Consequently, we can 
apply the Fokas-Gel'fand procedure with ODE (J5]) admitting a Lax representation (pi) 
given by the formula (35) with tangent vectors ^ for matrices A and S as in (16) and 

Two possible choices for a scalar product can be introduced on the tangents to the 
surface F e si(2,R) with the basis given by 



ei 



1 







1 




-1 


-1 


, e 2 = 


1 





, e 3 = 


1 



(48) 



In the first case, we decompose the matrix into the s£(2,R) basis and then use the 
standard Euclidean metric. The inner product and its norm in Euclidean space are 
given by 

\\X\\ = y/X'X*, (49) 



(X, Y) = X l Y\ 

where X = X l e u Y = Y^e j G sl(2, 
surfaces are Riemannian manifolds. 



hJ 



1,2,3. With the inner product (49) the 



In the second case, we use the Killing form on si (2, M), which is given by [3| [TT 

1 



B(X,Y) = -tr(XY). 



(50) 
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In terms of the basis (48), the matrices X, Y e sl(2,~. 
represented as follows 



B(X\Y j ) 



Bi 



1 
1 
0-1 



and the Killing form can be 

(51) 
(52) 



So, the Killing form has signature (2, 1) and induces a pseudo-Euclidean metric. The 
surfaces F 6 sl(2, R) are pseudo-Riemannian manifolds. 

Let us now explore certain geometric characteristics of the surfaces immersed in 
the s/(2,M) algebra. These geometric properties include the fundamental forms, mean 
and Gaussian curvatures. For example, the first fundamental form for the surfaces F 
and F Ux for any function f(u) with the pseudo-Euclidean metric are given by 



Ib{F st ) 



and 
Ib(F u * 



ff" 



f-9 f-sf 



u 



(u + X) 3 

Vf , V{f-g) 



dxdy+2 



u + X (u + X) 2 (u + Xf 



dxdy+ 



2 



v + X 



Vf 



f-9 

(u + A) 2 

VXf-g) 



u + x (u + xy 



dy\ (53) 



dy 2 .(5A) 



Note that in both cases, the first fundamental forms admit null vectors in the dx 
direction. 



4. The static (ft 4 field equation and its soliton surfaces 

We now present an example which illustrates the theoretical considerations. We intend 
to discuss in detail the construction of static and translation-invariant solutions of the 
4 -field equations [12] 

A B 

AM = — M + — M 3 , 0<AB,D,el, (55) 
2D 2D v ' 

where A denotes the Laplace operator on variables (x,y,z). For the purpose of this 

investigation, we limit ourselves to the three-dimensional Lie subalgebra spanned by 

{Li, P2, P3}, where the infinitesimal generators of rotation L\ and translations P2 and 

P 3 are 

d d d d 

L 1 = y --z-, P2 = ^, ^3 = ^- (56) 
az ay ay az 

The invariant solutions are of the form [T7] 

M(x) = u(0, e = e(x-f ), |ef = l, (57) 

where xo and e are constant vectors and the variable £ is obtained by applying the 
rotation L\ and translations P2, P3 to the symmetry variable £ = x. The translationally 
symmetric solution w(£) satisfies the second-order equation 

u xx = -2k 2 u 3 + (k 2 - h)u, (58) 
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-1 R -B 

h = —■ (59) 



2D v 2 ' 



AD 



Integrating once, (p8l) the first integral is 



(60) 



For specific choices of constants k\ and k 2 , we obtain different Jacobian elliptic functions 
sn cn and dn [U, |2] . 



fci fc 2 Solutions of ( 60 ) 



(61) 



1 —k 2 sn(x, k) 

k' 2 k 2 cn(x, k) 

—k' 2 1 dn(x, k) 

The moduli k of the elliptic functions are chosen in such a way that k' 2 + k 2 = 1 and 
< k, k' < 1. This ensures that the elliptic solutions possess one real and one purely 
imaginary period. 

The potential matrices L and M take values in the sZ(2,R) Lie algebra and are 
polynomial in u of third degree [5] 



M 



L 



u x (u-\)(k 2 (u 2 + \ 2 ) + k 1 -k 2 ) 
u + A — M x . 

-3£; 2 u 2 + 2AA; 2 m + k l -k 2 - k 2 X 2 

1 



G sZ(2,R) 
G sZ(2,R) 



(62) 
(63) 



where the discriminant is chosen as 

^(A) = /(-A) = (l-A 2 )(fc 1 + A; 2 A 2 ). (64) 

Solving the LSP Q with potential matrices given by (62) and (63), the most general 
form of the wave function is [9] 

(V / 9P0-^)*+-(v / 9W+« :c )*- (v / sW+« I )*--(a/9W-''x)* + 



$ = 

where 
ty± = exp 



2 v / m+A 

/u+A(^++^_) 
2 ' 



2 A /g(A)Vw+A 
yn+A(1'--^' + ) 

2V9(A) 



G SX(2, 



(65) 



f - uiu ^ 2 



+ c 



X 



A^i V ' \ 2 ' V h 

2^(X) v /(l-u 2 )(k 1 + k 2 u 2 ) + (k 2 -k 2 - 2k 2 X 2 )u 2 + (k 2 - h)X 2 + 2fci 



(66) 



_2 v ^(A)v / (l -M 2 )(fci + fc 2 u 2 ) - (A; a - fc 2 - 2A; 2 A> 2 - (fc a - fc x )A 2 - 2fc a _ 

where Co is a real integration constant and IT (u, a, b) is the normal elliptic integral of 
the third kind, see e.g. [2] 

dt 



II (u, a 2 , k) 



[l - a 2 t 2 )VT^¥VT^W 



(67) 
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In the graphs below, we choose the integration constant Co so that ^±(0, 0) = 1. 



The surfaces associated with the elliptic equations (58) are given by the formula 



(35). As an example, we consider separately three cases of the Jacobian elliptic function 



sn(x,k). When the discriminant g(X) > 0, the surfaces display exponential type 
behavior and when g(X) < the surfaces behave like trigonometric functions (see figure 

In the pseudo-Euclidean metric defined by the Killing form, the first and second 
fundamental forms for the surfaces F ST and with Q = u x are 

I(F ST ) = -k 2 {u - \)dxdy + (2k 2 u 2 - AXk 2 u + Qk 2 \ 2 - 2k 2 - 2) dy 2 
II{F ST ) = 2- 1 ' 2 k 2 {u - \)dx 2 - 2 1 ' 2 k 2 {u 2 - \ 2 )dxdy 

+2 3 ' 2 (jfeV - X k 2 u 2 + (k 2 X 2 -k 2 -l)u + k 2 X 3 ) dy 2 
I(F Q ) = k 2 (u 2 - l)(k 2 u 2 - l)(3u - X)dxdy 



2 ^ v 3 



+2 ^ fcV + 4 k 4 u 5 X - 2 k 4 u 4 X 2 - 4 k 2 X (l + k 2 ) 

+2 k 2 (-3 + A 2 + A; 2 A 2 ) u 2 + 4 A k 2 u + 2 - 2 A; 2 A 2 + 2 k 2 ^j dy 2 

= ^A: 2 (A -3n) ^ 2 + _ 3u)(u + x)dxdy 

v2 

+ 2V2(U + X) f 1 + + g ^ _ k 2_ 2Xk 2 u _ u 2 k 2 x 2 _ 2 fcVA 



+2 m 3 A k 2 - k 4 u 2 X 2 + 2 k 4 X u 3 + k 4 u 4 X 2 - k 4 u & ) dy 2 

The normals are 

(2k 2 u 2 ~k 2 ~l)uV2 



1 



N(F®) 



V2 
-1/2V2 

1/2^ 



and the Gaussian and mean curvatures are 

K(F ST ) = 2 (2 fcV - A; 2 - l) (A - u) A; 2 u 
#(F 5T ) = A; 2 A 2 - 2 A k 2 u + 3 fc V - k 2 - 1 

if(F Q ) = 2 A; 2 (ix + A) (A - 3 u) (-1 - k 2 - 3 k 2 u 4 + 6 fcV + 2 fcV - 3 fcV) 
H(F Q ) = y/2 (A - 3 it) u.A; 2 ^ - 5 fcV - 2 fcVA + A; 2 (6 + A; 2 A 2 + 6 A; 2 ) « 4 

+2 k 2 X (1 + A; 2 ) M 3 - A; 2 (9 + A; 2 A 2 + A 2 ) u 2 - 2 A A: 2 M + A; 2 A 2 + 1 + A^ . 

Note that the second fundamental form and the Gaussian curvature for the surface 
F ST in the pseudo-Euclidean metric coincide with those given for the surface in the 
Euclidean metric j9]. Also, the surface F lies in a plane in the moving frame defined 
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F 4 : A = 1.2, k = 0.5, g(X) < 



A = 0.5, k = 0.5, g(X) > 



Figure 1. Surfaces F 4 for u = sn(x, k) and a; and y € [- 
components of the immersion function in the basis ( 48 1 



The axes indicate the 



by conjugation with respect to the wave function $. Graphs of these surfaces can be 
found in [9]. 

Finally, we consider the surfaces associated with dilation symmetry (i.e. with 
constant in (35)). The surface is given by 

F i = (68) 

and the first fundamental form is 

I(F 4 ) = (3 /2 k 2 u 2 - k 2 (A - 3 xu x ) u + 1/2 k 2 \ 2 - u x xk 2 \ -I /2 k 2 -I /2) dx 2 
+k 2 x 2 (u l - 2)(k 2 u 2 - l)(3u - \)dxdy 

+2 x 2 (W + 2 k VA - fcVA 2 - 2 k 2 \ (1 + A; 2 ) M 3 

+k i (_ 3 + a 2 + A; 2 A 2 ) u 2 + 2 A k 2 u + 1 - PA 2 + ^ rfy 2 . 

The other geometric characteristics are directly computable but are too involved to write 
out in an illustrative fashion. Below we gives graphs of these surfaces in figures 1. 



5. Final remarks and future perspectives 

The main goal of this paper is to discuss an adaptation of the Fokas-Gel'fand procedure 
for constructing exact, analytic soliton surfaces associated with integrable ODEs 
admitting a Lax representation, as presented in [9j. In our investigation of the immersion 
formula for 2D-surfaces immersed in a Lie algebra we have proceeded in the following 
manner: 

1 We have shown, as in the PDE case [8], the problem of immersion requires the 
examination of conformal symmetries in the spectral parameter, gauge symmetries 
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of the LSP and generalized symmetries of the associated ODE model and its LSP 
(see the formula for immersion (12)). We have also demonstrated addition terms 
in the immersion formula, associated with translation in the independent variable 
x and dilation symmetry (see the extended formula for immersion (J35j) ) . 

2 We have constructed a Lax pair for a second-order ODE depending on an arbitrary 
function f(u) which includes, among others, the case of elliptic equations. 

3 We found explicitly the most general form of the wave function. 

4 From the general solution for the wave function, we have constructed 2D-surfaces 
in the Lie algebra s/(2,IR) by analytic methods for the general form of the ODE 
{u x ) 2 = f(u), with arbitrary f{u). 

5 Next, we identified a generalized symmetry of the considered ODE and showed that 
the immersion function can be explicitly integrated as in (35). 

6 The procedure was illustrated for translationally invariant solutions of the Landau- 
Ginzburg equation leading to the Jacobian elliptic equation. We have given some 
geometric characterizations such as the first and second fundamental forms of the 
surfaces as well as the mean and Gaussian curvatures. Additionally, we have 
provided graphs of the surfaces for different range of parameters leading to diverse 
types of surfaces. 

In the next stage of this research, using a group theoretical technique, the authors 
plan to consider soliton surfaces immersed in Lie algebras associated with the KdV and 
MKdV equations 

ut + u xxx + Qu x u = 0, it t + u xxx - Qu x u 2 = (69) 

respectively. The stationary states of these equations, when u t = 0, are given by 
Jacobian elliptic functions. Thus, a natural extension would be to compare the surfaces 
associated with PDEs as presented in j8] to the ODE surfaces associated with their 
stationary states (here the auxiliary variable y becomes t). For this purpose, we can use 
the exact solutions of the wave function (40) for stationary states to expand analytically 
for solutions in the neighborhood of stationary states and look for solutions of the wave 
function in the PDE case. 

Another avenue for investigation is the possibility of performing asymptotic analysis 
for the study of PDE surfaces using the ODE surfaces as approximations. A further 
investigation might be to consider how to apply the recurrence operator to generalized 
symmetries of the KdV and MKdV model in order to obtain recurrence relations for the 
surfaces. These and other issues will be addressed in our future work. 
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